//73- 3> SfS 



CASEFILE 

COPY 

ICSA 

INSTITUTE FOR COMPUTER SERVICES AND APPLICATIONS 

RICE UNIVERSITY 



I CSA-275-02 5-007 


A RANDOM NUMBER GENERATOR FOR 
CONTINUOUS RANDOM VARIABLES 
by 

Victor M. Guerra, R. A. Tapia 
James R. Thompson 
Dept, of Mathematical Sciences 
Rice University 


ABSTRACT 

A FORTRAN IV routine is given which may be used to generate random 
observations of a continuous real valued random variable. Jo use 
the procedure it is necessary to have as input only a table 

{ (x ,F(x )); i =0,1,2, ... k,k + 1 ) where x < x < . ..<x 

(i) (i) (0) (1) 

and F(x ) ~ 0 , F(x )~-'l . 
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A necessary part of any simulation study is the sampling of values 
of X , a random variable with cumulative distribution function F(X) . 

F may be specified either by a mathematical formula or empirically by 
a table of (X,F(X)) values. We consider here the special case where 
X is a real valued continuous random variable. Regardless of the 
specific form of F(X) , the probability distribution of the random 
variable Y = F(X) is the uniform distribution U[0,1] on the unit 
interval. 

There are a great many satisfactory pseudorandom number generators 
to obtain values from U[0,1] , e.g. [2] 

(1 ‘ 1) ■■ y n+l = 3141592 6 y .+ 27182818 (mod 2 35 ). 

The problem of sampling from F(X) may be looked at, then, as equivalent 
to finding an inversion formula 

(1.2) X = F" 1 (Y) . 

When it is not convenient to express F ^(Y) as a mathematical formula 
we may use an empirically obtained table (x , F(x^. ^,) ; 

x (l) < x (2) < •*■ < x (k)’ F(x (l) ) ~ 0 ’ F(X (k) ) ~ pluS a 8 °° d 
interpolation routine F^^(Y) to take random values from U[0,1] and 



obtain their approximate inverses as random samples from "nearly F(X)" . 

A. random number generator is advocated which uses an interpolation 
based on the highly stable piecewise cubic interpolation formula of 
Akima [1] . The user reads in a table [x^^, > 


< x.,. < ... < x._ . < x . The routine then computes and stores 

(0) (1) (k) (k+1 ) 

the coefficients of X = F ^ (Y) . Calls on the random number generator 
will first obtain a pseudorandom number y from U[0,lj . If 
y < F(x. . ) , the random number given is x. 1 . , if y > F(x., . , the 

(.1) (.-U v.k; 

random number given is x. . . If F(x n .) < y < F(x. . ) , the random 

(k; (.1; (,k; 

number given is F ^"(y) . A listing of the random number generator is 
given in Appendix I . 

Naturally, the quality of the random number generator proposed is 
determined by the ability of F^ to approximate a variety of cumulative 
distribution functions. In Appendix II , we show how the interpolation 
errors of F ^"(Y) compare with those of linear interpolation in the 
case of the standard normal, Cauchy and double exponential distributions 
with location 0 and scale 1 . Exact values of F ^(Y) are given for 


Y = .50(.02) ,98(.005) .995 . Interpolated values of F^ and linear 
interpolation are obtained at points between the mesh values and the 
relative errors computed and compared. 
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A PPE NDIX I 

; i RE L EA S E 1.1 M A I N 1 DA T E = 7 232 7 1.0 / 19/ A6 

C IN TER POL ATION O F THE INV ERSE NORML DISTR I BUT I ON 

COMMON N » P 1 (50) »P2( 50) » P 3 ( 50 ) *P4 { 50 ) »X(50) 

C N IS THE NUMBER OF INTERP OLATING KNOTS, N SHO ULD BE LESS OR E O TO 50 

N = 50 

CA LL IN I T 

C M IS THE NUMBER OF RANDOM NUMBERS TO BE GENERATED M SMALLER THAN 20 

M = 20 

CALL RNG { M ) 

CA LL EX I T . 

END ' 


TNI T 


date 


72327. 


G 1 REL EA SE 


10/19/416 


SUBROUTINE INIT . 

C ****************** 

C T HI S J5U B R O UTI NE GEN ERA TES THE PARAMETERS OF T HE CU BICS THAT I 
C INTERPOLATE THE DATA KNOTS 

C AR R^Y S P l * P 2 i.P3 .P A CO N TA l N THE P AR AME T E R S 

COMMON N , PI ( 50 ) ,P2( 50) ,P3( 50 ) , P4 < 50 j , X ( 50 ) 

100 FORMAT (2F12.6) 

200 FORMAT (2X.10F12.6) 

C READ DATA * # * A fc * # * £ # # # $ * £ >f< A A $ # $ * 

READ 100 , ( { X( I ) ,P4 ( I) ) . 1=1 ,N) 

C »* **»*»*»» REA D DAT A FINISHE S *» »*** »** »* 

N2 = N-2 

S1=(P4(2)— P4( 1 ) )/(X(2)— X(1 ) ) •- 

S2=( P4( 3 )-P4 < 2) )/<X(3)-X(2) ) 

S 3=( P4 ( 3 ) - P4 ( 4 ) ) / ( X ( 3 ) — X £4J 1 _____ 

DO 10 I=3,N2 

S4 = < P 4 ( I + 2 > - P41H 1) )/(X( 1+2) - X( D 1 ) ) 

W1 = A3SCS2-S1) 

W 3 = AOS ( S4 — S 3 ) 

IF (Wl*W3> 20,30.20 

20 . P3 ( I ) = ( W 3 * S 2 + W1 »S3 )/(W 1+W 3 ) 

GO TO 40 

30 P3(I) = «S*(S2 + S3) _ 

40 SI = S2 

S2 = S3 ; 

10 S3 = S4 

DO 60 I =3 , N 2 . 

A= XC 1 + 1 ) - X( I ) 

P_y.Il = J_P 3_{_I ) + P3( 1 + 1 ) - 2.0*( P4 ( I +- 1 ) - P4 ( J ) > /A) /(A *A ) 

P2( I )={ 3 . 0*( P4< 1+1 ) -P4 < I ) ) /A - 2.0*P3(li - P3(I+1))/A 

60 C ON TINUE 

RETURN 

: end ; 




1 _RE LEASE ! ,_1 


RNG 


D A T_E_ _=_7 2 3 2 7 1 0 /\ 9 /4 6 


SUBROUTINE RNG ( M ) . 

COMMON N .PI ( SO ) , P2( 50) ,P3( 50 ) ,P4 ( 50 ) , X( 50 ) 

R E AL * 8 PN (2 0 ) I 34 , I 35 ; 

100 FORMAT (I6.9F12.6) 

C * »»»»» RAN DOM NUMBER GE NE RATOR «* ****»» 

134 =57721566 

135 = 2 , 0**35 

00 1 0 1=1 . M 

I 3 4 = ( 31 415 926 *134 + 27182828 ) 

134 = ( DMOD( 134, 135) ) 

RN ( I ) = I 34/ I 35 ^ . 

PRINT 1 00 , I ,RN < I ) 

IF ( R N(I) » LE»X(3)) RN(1) = X(3) 

10 IF ( RN ( 1 ) . GE. X ( N-2 ) ) RN( I)= X(N-2) 

I_C _ =_ 1 

N3 = N-3 

DO 2 0 I =_3.N3_ ; 

00 30 J=1 . M 

V = R NJ J ) X ( I ) • 

IF ( R N ( J ) — X ( I ) ) 4 0.50.50 

5 0 IF (_R_N( J±_- _X_C_I + 1 ) ) 60.6 0,4 0 

60 YRN = Pt(I)*V*V*V + P 2 ( I ) * V *V + P3(I)*V + P4(I) 

C YRN I S T HE G EN ER ATE D RA NDOM OBSER VAT ION 

PRINT 100, IC ,RN( J) , YRN 

J IC = IC M 

40 CONTINUE 

I F { I C . GT . M ) GO TO 9 

30 CONTINUE 

_2_0 CO NT INUE , 

9 CONTINUE 

R ETURN 

END 




APPENDIX J I 


NOR 

M A L D I S 

T R I Q U T 

ION 

F ( X ) = 

X = 

E ( AK I M AS ) = 

E ( LINEAR )= 

0.50000 

- 0. 00000 


- -- — . 

0.50333 

0 . 00833 

0.00002 

0.00007 

0.50667 

0.01 666 

0.00003 

0.00007 

0.51000 

0 . 02500 

0. 00000 

0. 00000 

0.51 333 

0. 03334 

0. 00001 

0.0001 3 

0.51667 

0.04168 

0.00002 

0.00014 

0.52000 

0. 05002 

0. 00000 

0.000 0 0 

0.52333 

0 . 05337 

0.00002 

0.00020 

0.52667 

0 . 06672 

0.00002 

0.00021 

0.53000 

0.07508 

0. 00000 

0 . OCOOO 

0.53333 

0 . 00345 

0.00001 

0.00027 

0.53667 

0.09182 

0.00002 

0 . 00028 ' 

0.54000 

0 . 10020 

0.00000 

0.0 0 00 0 

0.54333 

0 . 10858 

0.00002 

0.00034 

0.54667 

0.1 1698 

0.00002 

0.00035 

0 . 5 5 0 0 0 

0 . 12538 

0.00000 

0.00000 

0 . 55333 

0 . 13379 

0.00002 

0.00041 

0.55667 

0. 14222 

0.00002 

0.00042 

0.56000 

0.15065 

0. 00000 

| 

O 

0 

O 

0.56333 

0. 15910 

0.00002 

0.00048 

0.56667 

0 . 16756 

0.00002 

0.00049 

0.57000 

0. 1 7603 

0. 00000 

0.0 

0.57333 

0 . 1845 1 

0.00002 

0.00055 

0.57667 

0 . 19301 

0.00002 

0.00056 

0.58000 

0.20152 

0. 00000 

6.6 

j 0.53333 

0.21 005 

0.00002 

0.00063 

0.58667 

0.21 859 

0.00002 

0.00064 

i 0.59000 

0.22715 

0.00000 

0.00000 

| 0.59333 

0 .23573 

0.00002 

0.00070 

0.59667 

0 . 24432 

0 . 00002 

0.00071 

i 0.60000 

0 .25293 

0 . oocoo 

0.00000 

0.60333 

0.26157 

0. 00002 

0.00078 

0 .60667 

0.27022 

0.00002 

0.00079 

0. 6 1 000 

0 .27889 

0.00000 

0.00000 

0 .61333 

0.28759 

0.00002 

0.00086 

0.61667 

0 .29630 

0.00002 

0.00087 

0.62000 

0.30504 

0.00000 

0.00 000 

0.62333 

0. 31 381 

0 . 00002 

0 .00094 

0.62667 

0.32260 

0 . 00002 

0.00095 

0.63000 

0.33141 

0.00000 

0.00000 

0.63333 

0. 34025 

0. 00002 

0.00103 

‘ 0.636 67 

0.34912 

0. 00002 

0.00104 

0.64000 

0 . 3580 1 

0. 00000 

0.00000 

0.64333 

0.36694 

0.00002 

0.001 1 2 

0.64667 

0.37589 

0.00003 

0.00113 

0 .65000 

0.38488 

0.00000 

0.0 

j 0.65333 

0.39389 

0.00002 

0.00121 

0 . 65667 

0 . 4 C 294 

0.00003 

0.00122 

< 0. 66000 

0.41202 

0. 00000 

0.0 

0 . 66333 

0.42114 

0.00002 

0.00131 

0.66 667 

0.43029 

0.00003 

0.00132 

0.67000 

0 . 43948 

0.00000 

0.00000 

0.67333 

0.44871 

0.00002 

0.00140 

0.67667 

0 . 45797 

0.00003 

0.00142 

0.68 0 00 : 

0. 46727 

0.00000 

0.00000 

0.68333 i 

0.47662 

0.00003 

0.00151 

0.68667 . 

0 . 4860 1 

0.00003 

0.00152 

0.69000 

0 . 49544 

0. 00000 

0.00000 ... 



■■ ' 


0.69333 

0.50491 

0. 00003 

0.00162 

0 . 69667 

0.51443 

0. 00003 

0.00163 

. 0.70000 

0 . 52400 

0. 00000 

0.00000 

0.70333 

0. 53362 

0.00003 

0.00173 

0 . 70667 

0.54328 

0.00004 

0.00175 

0 . 7 1000 

0 . 55300 

0.00000 

0.00000 

0.71333 

0.56277 

0.00003 

0.00185 

0. 71667 

0 .57260 

0.00004 

0.00187 

. . 0.72000 

0. 58248 

0.00000 

0.00000 

0 . 72 333 

0 . 59241 

0.00003 

0.00199 

0 . 72667 

0.60241 

0.00004 

0.00201 

0.73000 

0.61247 

0. 00000 

O 

0 

O 

0 .73333 

0.62259 

0.00004 

0.00212 

0.73667 

0 . 63277 

0.00005 

0.00214 

0.74000 

0 . 64302 

0. 00000 

0.00000 

0.74333 

0 . 65334 

0.00004 

0.00227 

0.74667 

0.66373 

0.00005 

0.00229 

0.76000 

0.67419 

0.00000 

0.0 

0.75333 

0.68472 

0.00004 

0.00243 

0 . 75667 

0 . 69534 

0 . 00006 

0.0024 6 

0.76000 

0 .70603 

0. 00000 

0.00000 

0.76333 

0. 71680 

0.00005 

0.00260 

0 . 76667 

0 .72765 

0.00006 

0.00263 

0.77000 

0 . 73860 

0.00000 

0.00000 

0.77333 

0.74963 

0.00005 

0.00279 

0 . 77667 

0.76075 

0.00007 

0.00282 

0.78000 

0.77197 

0.00000 

0.00000 

0.78333 

0.78329 

0.00005 

0.00299 

0.78667 

0 . 79470 

0. 00008 

0.00302 

0.79000 

0 . 80623 

0.00000 

0.00000 

0.79333 

0.81786 

0. 00006 

0.00321 

0 . 79667 

0.82960 

0.00009 

0.00324 

0.80000 

0.84146 

0 . 00000 

0 .00000 

0.80333 

0.85343 

0. 00007 

0.00344 

0.80667 

0 . 86553 

0.00010 

0.00349 

0.81000 

0 . 87776 

0.00000 

0.0 

0.81333 

0.890 1 2 

0.00007 

0.00371 

0.81667 

0.90262 

0.0001 1 

0.00375 

0.82000 

0.91 526 

0.00000 

0.0 

0.82333 

0.92805 

0.00008 

0.00400 

0 . 82667 

0. 94099 

0.00013 

0.00405 

0.83000 

0.954 1 0 

0.00000 

0.00000 

0.83333 

0. 96736 

0.00009 

0.00433 

0.83667 

0. 98080 

0.00015 

0.00439 

0.84000 

0.99442 

0. 00000 

0.00000 

0.84333 

1 .00823 

0 . 00009 

0.00471 

0.84667 

1 .02223 

0. 0001 6 

0.00476 

0.85000 

1 .03643 

0. 00000 

0.0 

0,85333 

1.05085 

0.00012 

0.00511 

0.85667 

1 . 06548 

0.00018 

0 . 0051 8 

0.86000 

1 . 08035 

0.00000 

0.00000 . 

0.86333 

1 . 09546 

0.00013 

0.00560 

0.86667 

1 . 1 1083 

0.00024 

0.00569 . . ; 

0.87000 

1.12646 

0.00000 

0.00000 • 

0.87333 

1 . 1 4237 

0.00014 

0.00616 

0 .87667 

1 . 1 5858 

0. 00028 

0.00626 

0.88000 

1.17509 

0.00000 

0.0 , 

' 0.88333 

1.19193 

0.00018 

0.00680 j 

0.88667 

1.20912 

0.00033 

0.00692 

0.89000 

1 . 22667 

0 . 00000 

0.00000 

0.89333 

1 .24460 

0.00019 

0.00759 _ i 



Et r L INEAR ) = 


CAUCHY DISTRIBUTION 

F(X) = X= E ( AK. I MAS ) = 

0*500 0 0 0.0 


0 .50333 

0.01047 

0 .00002. 

o.aooio 

0. d0667 

0.02095 

0.00002 

0.00012 

0.51000 

0.03 143 

o 

* 

o 

0.00000 ' 

0.51333 

0.04191 

0.00010 

0.00032 

0.5 1667 

0. 05241 

0. 00012 

0.00034 

U t jcUUU 

u*u o c y i 

*u*uouuu 

U#i>UUUU * ■ — 

0.52333 

0. 07344 

0.00005 

0 .00054 

0.5 2667 

0.08397 

0.00005 

0.00 056 - ••• ■ 

0.53000 

0.09453 

0 .00000 

0 .00000 

0. 53333 

0.1 0510 

- 0.00005 

0.00076 • 

0 . 53667 

0. 1 1570 

0. 0 0005 

0o0^078 
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0. 54333 

0. 13696 

0. 00005 

0.00098 

0 .54667- 

0. 14767 

0.00005 

0. 001 01 - - • 

0 . 55000 

0.1 5838 

0.00000 

0.00000 

— 0.55333 

0. 169 14 

0 . 0 0005 

0.00121 - - - 

0 .55667 

0.17993 

0.00006 

0.00123 
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— • -*■ 

0 o 56333 

0.20 164 

0. 0 0005 

0.00143 

0.56667 

0.21256 

-■ 0.00 006 

— 0.00146 

0.57000 

0.22353 

0.00000 

0.00000 

0.57333 

0.23455 

0.00006 

0.00166 - - - 

0.57667 

0.24562 

0 .00006 

0.00169 

U o joUUU 

Uo cuO f D 

U • U U L)UU 

U • UlJJu U U 

0 o 58333 

0.26795 

0.00005 

0.00190 

. o . 5 8667 

0.27921 

0.00006 

0.00192 ' 

0.59000 

0. 29053 

0.00000 

0.0 

0.59333 

0.30192 

0. 00006 

0.002 13- — - - 

0.59667 

0.31338 

0 .00007 

0.00216 

0 boOOOO 

0 ♦ 

0 • 0 0 0 0 0 

0 o 0QG 0 0 

0.60333 

0.33654 

0.00006 

0. 00237 

0.6 0667 

0 #34824- 

0*00 00 7 

- 0 . 00240 

0.61000 

0.36002 

0. 00000 

0.0 

0.61 333 

0 • 3 7 190 

0.00006 

0 . 00262 

0.6 1667 

0. 38386 

0.00008 

0.00266 

0 ®o2000 

0.39593 

0.0 000 0' 

0*0 

0 .62333 

0 .40809 

0.00006 

0.00287 

t - - 0.6266 7 

0.42036 

— “ 0.0 0008 

- 0.00291 - 

0 .63000 

0.43274 

0. 00000 

0. 00000 

0 . 63333 

0.44523 

0 .00 00 Z' 

0 . 0031 3 

0.63667 

0.45784 

0.00009 

0.00318 

0 •640 00 

0.47056 

0*00000 

' O'. 0000 0 ' 

0.64333 

0.48342 

0.00007 

0.0034 0 

0.64667 

0.49640 

0.00010 

0.00344 ' 

0 .65000 

0.50953 

0.00000 

0.00000 

0.65333 

0.52279 

0.00008 

0 .00367 ' - 

0 • 65667 

0. 53620 

0.00010 

0. 0 0372 

0 .6 60 00 

0.54975 

0 . U 0 oou 

0*0 - 

0.66333 

0. 56347 

0.00008 

0.00396 

0 . 66667 

0.57735 

0.00011 

0.00400 ' — ' ' 

0.67000 

0.59140 

0.00000 

0.00000 

0.67333 

0.6 0562 

0.00009 

0.00425 • 

0 .6 7667 

0.62003 

0.00012 

0. 00430 

0.68000 

0 .63462 

~ 0.0000 0' 

0.0 

0.68333 

0.6494 1 

0.00009 

0.00456 

0.68667 

" 0.66440 

0.00013 

0.00462 ' — ' 

0.69000 

0.67960 

0.00000 

0.0 


0.69333 

0.69502 

0 .0001 0 

&* 0 04 87 

0.6966 7 

0.7 1066 

0.00014 

0*O'O493 ~ 

0. 70000 

0. 72654 

0.00000 

0.0000 

0 .70333 

0.74267 

0 • 0 0 0 1 i 

0*00521 

0.70667 

0.75904 

0.00015 

0*00527 

_ 0#7 iooo 

0. 77568 

o.ooooo 

0*00000 — 

0 .7 1333 

0. 79259 

0.00012 

0.0 0556 

0.71667 

0. 80978 

0.0 0017 — - 

0«0 0563 

0.72000 

0. 82727 

0.00000 

0*0 00 0 0 

0 • 7 2. 3 3 3 

0.84507 

0 • 0 0 0 l 3 

0*00592 ' 

0. 72667 

0.863 1 8 

0.00019 

0*0 0600 

~ 0.73000 

0.88162 

— 0. 00000 " - 

Oo» 0 

0.73333 

0.90040 

0.00014 

0*0 06 31 

‘ 0. 73667 

0. 91 955 

- 0.00020 

0*00640 

0 .74000 

0.93906 

0.00000 

\0m 0 00 00 

0 • 74333 

“0T9 5 8 97" 

O .00015 

i© o*0 uo / 2 

0 . 74667 

0.97927 

0.00022 

0*0068 1 

0.75000 

liOOOOO 

- " 0.00000 

0*00000 

0.75333 

1.02117 

0.00015 

0*0 0 716 

0 . 756 6 7 

1. 04279 

“ ----- 00 025 — 

0*00726 

0.76000 

1 .06489 

0.00000 

:0»O 0000 

0.7 6333 

1 • 08749 

0*00017 

0 a# 0 0 7 6 1 

0.7 6667 

1.11061 

0.00029 

10*0 07 74 

~ 0.77000 

1.13428 

' 0 .0000 0“ 

Q* 0 

0.77333 

1.15851 

0.00019 

0*0081 1 

0 . 77667 

1 . 1 8 3 34 

0.00 03 1 

0*0 08 23 

0.78000 

1 .20879 

0 .00000 

0*0 0000 

Oo 7 8333 

1 V 2 3 4 9 0 

0.00021 •" “ 

0*0 0 865 

i 0 .78667 

1.26169 

0.00034 

0*00878 

, '• 0.79000 

l .28919' 

0.00000 

0*0 ' 

0.79333 

1.31745 

0.00023 

0*00922 

• 0 .7 96 67 ' — 

1 .34650 

0.00039 

0-00937 ' 

• 0.80000 

1 .37638 

0.00000 

0.0 

o .30333 

1.40714 

0 « 0 0 025 

0 to 0 0 ^ 8 3 

0.80667 

1.43881 

0.00044 

0*0 1001 

0.81000 

1. 47146 

o.ooooo 

0*00000 ■ 

0 .3 1333 

1.50512 

0. 00028 

0.01053 

O.Q 166 7 

1.53987 

0.0 0 049 

© m 0 1 0 7 1 

0.82000 

1.57575 

0.00000 

© 4P 0 

0 © 8 23 33 
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0.0 0 030 

0*0 1128 i 

0.82667 

1.65120 

0.00056 

0*01148 

0.83000 
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